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ABSTRACT 

C*") ■ Semi-analytical models of galaxy formation can be used to predict the evolution of 

the number density of early-type galaxies as a function of the circular velocity at the 
virial radius, w c ,vir- Gravitational lensing probability and separation distribution on 
the other hand are sensitive to the velocity dispersion (or circular velocity) at about 
the effective radius. We adopt the singular isothermal ellipsoid (SIE) lens model to 
estimate the velocity dispersion at the effective radius. The velocity dispersion from 
' strong lensing based on the SIE, asm is then closely related to the observational 

, central stellar velocity dispersion, cr ccn t; we have empirically asm ~ c ccn t. We use 

radio lenses from the Cosmic Lens All-Sky Survey and the PMN-NVSS Extragalactic 
Lens Survey to study how the velocity dispersions, asm, are related to v c y lr ; if the 

galaxy were a singular isothermal sphere up to the virial radius, v c ,vir = V^csie- When 
we include both the lensing probability and separation distribution as our lensing 

constraints, we find ct S ie/(200 km s" 1 ) = [(1.17l^)w c , v ir/(200 km s" 1 )] ' 22 ^ °ol for 
200 km s -1 < asm ^ 260 km s _1 ; at ctsie = 200km s -1 , the ratio v / 2o"SiE/«c,vir is 
about 1.65lg;|? (68% CL) but decreases to 0.65±^ (68 % CL) for a S m = 260 km s -1 . 
These results are consistent with those of ISeliakl (|2002t ) obtained from galaxy-galaxy 
^ ■ weak lensing for galaxies of around L*. However, our results clearly suggest that the 

ratio must vary significantly as asm is varied and are marginally discrepant with 
the Seljak results at ctsie = 260km s^ 1 . The scaling asm "cv? 005 i s broadly 
consistent with those from galaxy occupation statistics studies and the most recent 
galaxy-galaxy weak lensing study. We discuss briefly the implications of our results for 
galaxy formations and structures. These constraints can be significantly strengthened 
when larger lens samples become available and the accuracy of semi-analytical model 
predictions improves. 

Key words: gravitational lensing - galaxies: formation - galaxies: haloes - galaxies: 
kinematics and dynamics - galaxies: structure 



1 INTRODUCTION 

The ACDM structure formation model has been very 
successful in explaining many observations, for example 



the cosmic microwave background radiation dSpergel et al 



2003 
2001 



2006) , and the large -scale structures l|Peacock et al 



Tegmark et al.l [20041 ). In this model, the tiny quan- 



tum fluctuations in the early universe are amplified due 
to gravitational instabilities, and eventually evolve into 



* chae@sejong.ac.kr 



highly nonlinear structures. The dissipationless cold dark 
matter dominates the gravity, and determines the ba- 
sic abundance, the internal structure and formation his- 
tory of the nonlinear dark matter haloes. Such infor- 
mation can be reliably obtained using JV-body simula- 
tions and/or an alytical models such as the Press-Schechter 
form alism (e.g., iPress fc Schechterl fl97i ISheth fc Tormenl 
2002). For example, the dark matter haloes follow ap- 
prox imately the Navarro , Fren k and White (NFW) pro- 
file l|Navarro et al.l fl997l , 12004 ) . and their sh apes can be 
approximated by tri-axial ellipsoidal models l|jing fc Sutol 
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2002). Their time evolution can be studied numerically, 
which is also matched well by the analytical extended Press- 
Schechter formalism (e.g.. | Bond et al.l Il99ll : lLacev fc Cold 
ll994l ; ISheth fc Tormenll2002l ). 

On the other hand, the visible galaxies we see today are 
formed by the condensation of baryons within dark matter 
haloes due to radiative processes. When the baryons sink to- 
ward the centre of dark haloes and become self-gravitating, 
star formation, active galactic nuclei and feedback pro- 
cesses will occur. Unfortunately, these phenomena are not 
well-understood. For clusters of galaxies, baryonic cooling 
is not very important, so their mass profiles still approxi- 
mately follow the predicti ons from iV-body simulati ons, as 
seen fro m kinematic (e.g.. van der Marel et"al1l2000h. lens- 
ing (e .g. JComerford et al.1l2006l ) and X-ray ( Voigt fc Fabianl 
2006) studies of clusters of galaxies. On the scale of galax- 
ies, cooling is important, and the central profiles of galaxies 
are likely to be significantly affected by baryonic processes. 
Hydrodynamical high-resolution simulations can now sim- 
ulate regions of the universe or individual galaxies (e.g., 
iMeza et alJl2003h . but they cannot yet resolve the inter- 
nal structure of galaxies and at the same time simulate 
large enough volume to be statistically representative (and 
realistic). Nevertheless, rapid progress has been made, in- 
cluding more reali stic modelling of the mul ti-phase in- 
terstellar medium (jSpringel fc Hernauistl 120031 ) . An alter- 
native approach is to use semi-analytical studies which 
can i ncorporate many physi c al processes in an intuitive 
wav teauifmann et all Il999l: [Somerville fc Primackl Il999l : 



ICole et all l2000l : iKang et all 120051 ; ICroton et all l2006h 
These models can reproduce the luminosity function and 
correlation function of galaxies reasonably well. However, 
although these models can provide the number of galaxies 
as a function of circular velocity at the virial radius, they 
cannot yet reliably predict the central dynamical properties 
of galaxies, such as the central velocity dispersions. 

The central velocity dispersions (or the velocity disper- 
sions in the optical regions) of galaxies are precisely the 
information needed for the studies of gravitational lens- 
ing, as the lensing probability is proportional to a 4 while 
the i mage separation is proportional to tr 2 . E arlier studies 
(e.g., iNaravan fe Whitelll98Sl ; iKochanekl [l995h used singu- 
lar isothermal spheres to constrain variants of CDM models. 
In this scenario, the velocity dispersion is simply related to 
the circular velocity at the virial radius by v c = \/2a. How- 
ever, this assumption is likely to be invalid for two reasons. 
The NFW profile for haloes does not give rise to a flat ro- 
tation curve. In fact, the value of its peak is about 20 per 
cent higher than that at the virial radius for a concentration 
parameter (~ 10) appropriate for Milky- Way sized haloes. 
Second, the baryons settled at the centre may further in- 
crease the central rotations and velocity dispersions. As the 
assembly of galaxies is not yet well understood, how the ve- 
locity dispersion in the optical region relates to the virial 
circular velocity is therefore a parameter we want to extract 
from observations. In this paper, we propose to use strong 
gravitational lensing to constrain this key parameter. 

At the time of this writing no study has made use of 
the predictions of semi-analytical studies of galaxy abun- 
dances of different types. Several studies have used gravi- 
tational lensing to constrain the evolution of galaxies dMao 



120031) . In this paper, we will use the predicted abundance 
of ellipticals from semi-analytical models as a function of 
virial circular velocity. We then adopt a simple form to pa- 
rameterise the relation between the velocity dispersion at 
the optical radius and the circular velocity at the virial 
radius. We use radi o lenses from the Cosmic Lens All - 
Sky Survey (CLASS: iMvers et al.|[2003l ; iBrowne et alj|2003l ) 
and the PMN-NV SS Extragalactic Lens Survey (PANELS; 
IWinn et al.ll2001bl ) to constrain this factor empirically. This 
in turn will provide an important way of understanding the 
baryonic effects on the central properties of galaxies. Our 
approach is independent of the galaxy-ga laxy weak lensing 
studies of the same factor |Scliak 2002]), who found that 
for both early-type and late-type galaxies around L*, their 
peak velocities are about a factor from 1.7 to 1.8 (with an 
uncertainty of about 20%) of the circular velocity at the 
virial radius. The structure of this paper is as follows. In §2, 
we outline our method and the data to be used. In §3, we 
present the results of our analyses. Finally, in §4 we discuss 
the implications of our results for galaxy formation along 
with possible sources of systematic errors and give the main 
conclusions of our study. 



2 MODEL, METHOD AND DATA 

In §2.1 we introduce velocity functions of galaxies start ing 
from the Schechter luminosity function jSchechtcr 1973) of 
galaxies. We also adopt a model for the relation between 
the velocity dispersion arising from strong lensing and the 
virial circular velocity predicted from semi-analytical mod- 
elling of galaxy formation (§2.1). In §2.2 we briefly describe 
our adopted semi-analytical model (SAM) of galaxy forma- 
tion and present simulated data to be used in this work. 
We also outline our formalism of strong lensing statistics to 
constrain the radial velocity profile of galaxies based on the 
SAM circular velocity function and strong lensing data. 



2.1 Velocity functions of galaxies 

The number density of galaxies as a function of luminosity 
is described by the Schechter luminosity function (LF) </>l 
given by 



dn = ^{L^dL = 4>* ( y— 



dL 



(i) 



We assume an effective power-law relation between the lu- 
minosity (L) and the velocity dispersion (a) for the scale 
under consideration of galaxies given bjQ 



L 
1Z 



/5 V D 



(2) 



Il99ll ; iMao fc Kochaneklll994l : lOfek et al.ll2003l:lchae fc Mao 



Then the number density of galaxies as a function of veloc- 
ity dispersion can be described by the velocity dispersion 
function (VDF) i^vd given by 



In reality the relation between L and a does not follow an exact 
power-law relation but suffers from a significant scatter. Thus, the 
power-law fitting function between L an d a, namely the Fabcr- 
Jackson relation l|Faber fc Jacksonlll976h . may not be the same 
as this effective relation (see lSheth et al] |20031. 
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Table 1. Summary of velocity and related parameters 



the virial circular velocity function (CVF) 0cv takes the 
following fornQ 



parameter meaning 



relevant scale 



v c v j r circular velocity at the virial radius r v ; r 

v c e fi circular velocity at the effective radius r e ff 

CT mot j lens model velocity dispersion ~ r c g 

(""sib) (SIE model velocity dispersion) 

Tccnt observed central velocity dispersion < 

V^tsie SIE model circular velocity 

V = CT mo d /(Tccnt 

/ = CT mod/^c,vir 



Notes. r v 



virial radius. r e ff = effective radius. 



dn — 0vD(o")dcr = 0» I — 



exp 



/3VD 



da 

PVD 

cr* 



where q V d = (o?l + 1)/3vd- 

In this work we shall study how the velocity disper- 
sion in the optical region and the inner halo constrained 
by strong lensing statistics is related to the circular veloc- 
ity at the virial radius (f c ,vir) predicted by an up-to-date 
semi-analytical model of galaxy formation. Notice that the 
velocity dispersion constrained by strong le nsing depends 
on the adopted lens model (e.g., IChaell2003l ). We shall de- 
note the lens model-dependent velocity dispersion by (j mo d. 
Parameter a moc j corresponds to the velocity dispersion ap- 
proximately at the scale of the optical region and the inner 
part of the halo because the observed galactic-scale lens im- 
ages are formed at around the effective radii. In particular, 
although a lens model may imply a mass distribution from 
the galactic centre to infinity, the mass distribution well out- 
side the optical region is an extrapolation that is not sensi- 
tive to strong lensing. We expect that the velocity dispersion 
o"mod is closely related to the central stellar velocity disper- 
sion, er ccnt that is observable in the aperture-limited spectro- 
scopic observations. In this work we shall assume a constant 
proportionality between the two parameters; we write 



(4) 



(see §4 for a discussion on the possibilities of varying rj). 

Hereafter for the sake of simplicity parameters a and 
v c shall refer respectively to a moc j and w c ,vjr unless specified 
otherwise. We parameterise the relationship between the two 
parameters a and v c by 



= f( v c)Vc 



(5) 



This is reasonable as on the cluster scale, the baryonic ef- 
fect may be small, while the baryonic condensation within 
galactic-sized halos may significantly boost the velocity dis- 
persion in the optical region relative to the virial circular 
velocity. We adopt a model given by 



f(Vo) = f. ^ 



(6) 



~dn — 4>cv(v c )dv c = 4>* 



«cv- 1 



exp 



/3cv 



where we have the following relations: 

acv = qvd(m+1)i Pcv = /3vd(m+1)> and cr» = /*w c *-(8) 

The above CVF (equation [7]) implies an effective power-law 
relation between the luminosity and the virial circular ve- 
locity of galaxies given by 

/3 CV 

(9) 



L_ 

17 



The velocity and related parameters introduced above are 
summarised in Table 1. 
(3) Finally, a galaxy of v c * may not correspond to a typical 
L» galaxy. Suppose an L* galaxy has a velocity dispersion 
of cto- Then the value of / (equation O for ctq is given by 



/o = /» 



/V(H-l) 



(10) 



We shall use parameter fo rather than /* for a chosen value 
of (TO. 



2.2 Strong lensing statistics and semi-analytical 
model of galaxy formation 

The statistical properties of strong lensing are determined 
by four ingredients, namely the number density of potential 
lenses as a function of redshift, the lens cross sections, the 
cosmology, and the magnification bias depend ing on the po- 
tential source population (see, e.g.. lChaell2003l ). Most studies 
of strong lensing statistics have used the present-day galaxy 
population given in the form of luminosity (or velocity dis- 
persion) function for the number density of potential lenses 
(sec, e.g., Turner et al. 1984; Fukugita ct al. 1992; Kochanck 
ll996tlHelbig et al.lll999l : IChaell200j ; lMitchell et al.ll2005l ). In 
this paper, however, we adopt the numbe r density of galaxie s 
from the recent semi-analytical model of iKang et al. {2005 ), 
which gives the number density of galaxies as a function of 
redshift and the circular velocity at the virial radius. 

The galaxy catalogue constructed bv lKang" et al.l (l2005f) 
used a high-resolution N-body simulation by Ijing fc Sutol 
(2002). The simulation follows the evolution of 512 3 particles 
in a cosmological box of 100/i~ 1 Mpc with the standard con- 
cordance cos mological model (fl m ,n = 0.3, CIa,o = 0.7, as = 
0.9, h = 0.7). Hang etafl (|2005h adopted reasonable treat- 
ment of physical processes in predicting galaxy population 
such as cooling, star formation rate and energy feedbacks 
from supernova, and dust extinction. The parameters in the 
model were normalized by the local galaxy luminosity func- 
tions. Th e readers are refere ed to that paper for more detail. 

The IKang et al. I (120051 ) SAM reproduces the luminos- 
ity function of galaxies reasonably well for the range of 
luminosity suitable for strong lensing studies. The semi- 
analytical studies provide a simple way of dividing galaxies 



where n — corresponds to the case where the velocity 
dispersion is related to the virial circular velocity by a con- 
stant proportionality. For the above model of equation ((6| 



2 This particular form is the same functional form as the VDF 
given by equation J3). This is the direct consequence of adopting 
the particular model of equation IBT l. 
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Tabl e 2. Distrib ution of the number of early-type and late-type galaxies as a function of virial circular velocity predicted by the 
Kang et al ] l|2005h semi-analytical model of galaxy-formation. Each value given here is the number of early-type (late-type) galaxies in 
a velocity bin of width 25 km s — 1 contained within a total comoving volume of 10 X 10 6 h~ 3 Mpc 3 distributed for the redshift range of 
0.3 <Z < 1. 



% 

(bin #) 


f c ,vir(«) (km s 1 ) 

(mean virial circular velocity) 


AA r(carly) (j) 

(# of early-type galaxies) 


AA T(late)(j) 

(# of late-type galaxies) 


1 


112.5 


30321 


148585 


2 


137.5 


12854 


73322 


3 


162.5 


5816 


40381 


4 


187.5 


2789 


22903 


5 


212.5 


1719 


14661 


6 


237.5 


1377 


8802 


7 


262.5 


1266 


5899 


8 


287.5 


830 


3808 


9 


312.5 


763 


2357 


10 


337.5 


688 


1484 


11 


375. 


532 


830 


12 


425. 


295 


442 


13 


475. 


186 


260 


14 


525. 


145 


120 


15 


575. 


108 


72 



Table 3. Strongly lensed systems f rom the CL ASS jBrowne et alj|2003f) and the PANELS southern sky jWinn et al.ll2000l . l2001al lbl. 
l2002al lbh radio surveys (taken from Chac 2005[ ). The systems under the 'CLA SS statistical' survey are members of the well-defined 
CLASS statistical sample of 8958 radio sources <Browne et al.ll2003l ; IChaell2003l) . 



Source Survey Source Lens Image Separation Image Lens Type 

Redshift Redshift (arcsec) Multiplicity 



B0128+437 


CLASS 


3.124 


1.145 


0.54 


4 


unknown 


J0134-0931 


PANELS 


2.225 


0.7645 


0.681 


2+4 


2Gs 


B0218+357 


CLASS statistical 


0.96 


0.68 


0.334 


2 


spiral 


B0414+054 


CLASS 


2.62 


0.958 


2.09 


1 


early-type 


B0445+123 


CLASS statistical 




0.558 


1.33 


2 


early-type 


B0631+519 


CLASS statistical 




0.620 


1.16 


2 


early-type 


B0712+472 


CLASS statistical 


1.34 


0.41 


1.27 


4 


early-type 


B0739+366 


CLASS 






0.54 


2 


unknown 


B0850+054 


CLASS statistical 




0.588 


0.68 


2 


spiral 


B1030+074 


CLASS 


1.535 


0.599 


1.56 


2 


early-type 


B1127+385 


CLASS 






0.70 


2 


spiral 


B1152+199 


CLASS statistical 


1.019 


0.439 


1.56 


2 


unknown 


B1359+154 


CLASS statistical 


3.235 




1.65 


6 


3Gs 


B1422+231 


CLASS statistical 


3.62 


0.34 


1.28 


4 


early-type 


B1555+375 


CLASS 






0.42 


4 


unknown 


B1600+434 


CLASS 


1.57 


0.415 


1.39 


2 


spiral 


B1608+656 


CLASS statistical 


1.39 


0.64 


2.08 


4 


2Gs 


J1632-0033 


PANELS 


3.42 


1 


1.47 


2 


early-type 


J1838-3427 


PANELS 


2.78 


0.36 


1.0 


2 


early-type 


B1933+503 


CLASS statistical 


2.62 


0.755 


1.17 


4 


early-type 


B1938+666 


CLASS 


> 1.8 


0.881 


0.93 


ring 


early-type 


J2004-1349 


PANELS 






1.13 


2 


spiral 


B2045+265 


CLASS statistical 




0.867 


1.86 


4 


puzzling 


B2108+213 


CLASS 




0.365 


4.55 


2 or 3 


2Gs+clustcr 


B2114+022 


CLASS statistical 




0.32/0.59 


2.57 


2 


2Gs 


B2319+051 


CLASS statistical 




0.624/0.588 


1.36 


2 


early-type 



into two populations, namely the early-type (ellipticals and 
SO's) population and the late-type population based on the 
correlation between the B- band bulge-to-disk ratio and th e 
Hubble type of the galaxy (|Simien fc de Vaucouleurs| [l986). 
We use only the early-type population for this work because 
the radio lens sample we shall use contains only a small num- 



ber of spiral lensing galaxies. Table 2 shows the numbers of 
early-type and late-type galaxies in velocity bins contained 
in 10 x 10 6 h~ 3 Mpc 3 volume distributed for the redshift 
range of 0.3 < z < 1. 

While the SAM data can give a virial circular velocity 
function, strong lensing is not sensitive to the virial circu- 
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lar velocity but the velocity dispersion in the optical region 
and the inner halo. Hence to proceed in lensing computa- 
tions using the number density of early-type galaxies from 
the SAM we must relate the virial circular velocity to the 
velocity dispersion. We proceed as follows. We adopt the 
VDF (^vd(ct) (equation [3} for lensing computations and the 
CVF 4>cv(v c ) (equation [7]) for describing the SAM data as- 
suming that they are related by the power-law model given 
by equation ([6jl . We then fit lensing data and the SAM data 
simultaneously. To fit the VDF to lensing data we use a 
ma ximum likel i hood method based on t hose r e cently used 
by IChae et al l l|2002t ). IChad [|2003h and IChad l|2005l ). The 
likelihood for lensing is defined by 

in c = rr wi in 5 PXS (j) + Mr mi - p(*o] + E ln mo .(h) 

\j=1 / \k=l 1 = 1 J 

where the first term is the likelihood due to the relative 
im age separatio n probabilities of lensed sources as defined 
by IChad (|2005h and the terms in the second parenthe- 
sis are the likelihood due to the C LASS statist i cal sa mple 
dBrowne et al.l l2003h as defined by IChae et all (|2002h and 
IChad [|2003h ~ p arameters W k are the weight factors, 5pm(j) 
(equation 4 of IChad I200I ) are the relative image separa- 
tion probabilities for lensed sources with well-defined image 
separations (Ta ble 3; see a lso Table 1 of IChadl2005l ). p(k) 
(equation 40 of IChaell2003h are the total lensing probabili- 
ties for the unlensed s ources in th e CLASS statistical sample 
(see §3.1 a rid 53.2 of IChadl2003l ). and Sp(l) (equations 29, 
38, or 39 o f lChadl2003¥ are the differential lensing probabil- 
ities for the lensed sources in the CLASS statistical sample 
(Table 3; see also Table 1 of I Chad 120031 '). Notice that the 
above lensing probabilities are calculated assuming the sin- 
gular isothermal ellipsoid ( SIE) lens m odel whose projected 
surface density is given by l|Chadl2003h 



£sie(:t, y) = 



°SIE 



v^A(e) 



2G x 2 + e 2 y 2 



(12) 



where ctsie i s the model velocity dispersion, e is the axis ra- 
ti o, and A(e) i s the 'dynamical normalisation' factor (§ 2.1.1 
of IChad 12003 ). We assume that galaxies are not biased to- 
ward either oblate or prolate shape and have a mean pro- 
jected ellipticity of e(= 1 — e) < 0.45 so that A(e) « 1. Notice 
that the shape and characteristic velocity dispersion of the 
VDF are most directly constrained by the image separation 
distribution. The absolute lensing probability of the CLASS 
statistical sample is also sensitive to the VDF as the lens- 
ing probability is proportional to 0,cr*. The chi-squared for 
lensing is then defined by 



Xlc 



-2m£. 



(13) 



To fit the CVF (equation [7} to the SAM data (Table 2) we 
define the following chi-squared function 



2 ( <t>cv(.v c (i))Av c AV - AN(i)^ 2 

XSAM — 2_^i 



(14) 



where Av c = 25km s^ 1 and AV = 10 7 h' 3 Mpc 3 . Finally, 
the total chi-squared function is defined to be 



222 

Xtot = Xlens + XSAM ■ 



(15) 
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Figure 1. The numbers of early- type galaxies and late-type 
galaxies (Table 2) as a function of circula r velocity at the virial 
radius predicted by the iKane et al. 1 1120051 ) semi-analytical model 
of galaxy formation. The circular velocity functions fitted to the 
data using equation J7J are also shown. We have used Poisson 
errors for the fitting, which are however too small to display, and 
ignored the first four data points of the early-type population 
(see the texts). The results are for 10 X 10 6 h — 3 Mpc 3 comoving 
volume nearly uniformly distributed over the redshift range of 
0.3 < z < 1. The velocity bin size is 25 km s _1 . 



3 RESULTS 



We fir st examine the simulated data from the iKang et al.l 
|2005t ) SAM. Fig. [Q shows the SAM data points (Table 2) 
along with the fitted CVFs (equation [TJl through the chi- 
squared function given by equation (|14p. Notice that the 
lowest 4 data points for the early-type population are not 
included in the fitting because they do significantly deviate 
from the behaviour of the rest of the data points. We shall 
not use these 4 data points in our study of the velocity profile 
of early- type galaxies (see below and §4). Fig. [2] shows the 
constraints on the shapes of the CVFs and compares with 
th e shap e of th e VDF of the early-type population obtained 
bv lChad l|2005h . Notice that the CVF shape is significantly 
different from the VDF shape for the early-type population. 
This means that the velocity dispersion cannot be related 
to the virial circular velocity by a constant proportionality 
(i.e. the case /j, = in equation [6jl . 

We simultaneously fit the VDF (equation [3jl and the 
CVF (equation [7]) for the early-type galaxy population as- 
suming the relation between the two given by equations <(Sj 
and ([6]) to the lensing data and the SAM data by min- 
imising the total chi-squared function Xtot given by equa- 
tion (|15|) . We do not use the first 4 data points (i = 1 to 
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Figure 2. Constraints on the velocity functions of galaxies. The red and blue contours are resp ectively the c onfide nce limits (CLs) 



on the CVFs (equation for the early-type and the late-type populations as obtained from the Kan g et al.l 
data. The black contours are the CLs on the VDF (equation [3j for the early-type population as obtained by 



SAM simulated 
2005h from the 

image separation distribution of gravitational lenses. The three contours on each function correspond respectively to la, 2cr, and 3<r. The 



filled triangle represents the SDSS VDF as obtained bv lSheth et al.l d2003T) f r om th e measure d central veloc i ty dis persions of early-type 
galaxies. The filled square represents the SSRS2 VDF as obtained by IChaj l|2003h from the lMarzke et ahl l|l99Sl ) early-type LF and a 
measured Faber-Jackson relation. 



4 in Table 2) for the early-type population from the SAM. 
These 4 data points do not fit into the model CVF (equa- 
tion [7|) while all the rest do (Fig. [IJ. These data points 
m ight be an artifact due to the currently imperfect SAM 
of iKang et all (|2005h (see §4). Perhaps, more importantly, 
lowest velocity galaxies may not be effective for strong lens- 
ing for the following two reasons: the lensing cross-section 
scales as a STF , and their inner mass pro files may be too sha l- 
low (e.g. Jde Blokll2005l . see, however, lijwaters et alj|2003h . 
So simply ignoring them would make small errors in our 
study. Nevertheless, we must bear in mind that any of our 
derived results for « c ,vir < 190 km s^ 1 will be an extrap- 
olation and the strict range of validity for our analyses is 
200km s" 1 < « c , vir < 600km s" 1 (see Table 2). 

We consider first the case where we use both of the 
parts in equation {TJJ, namely both the relative image 
separation probabilities and the lensing properties of the 
CL ASS statistical sample. For the CLASS statistical sam- 
ple (|Browne et al.ll2003l ; IChaell2003l ) the number of sources 
that are not strongly lensed is Nv = 8945. There are 
13 strongly lensed sources in the CLASS statistical sam- 
ple as can be found in Table 3. However, the lenses for 
B0218+357 and B0850+054 are spiral galaxies and the sys- 
tem B2045+265 might include a spiral galaxy (or spiral 
galaxies) but its len sing interpretation remains puzzling (see 
iFassnacht al.|[l999l ). Excluding these three sources we take 
JVl = 10 for lensing solely due to the early-type galaxy 
population. To be included in the likelihood as relative im- 
age separation probabilities are B0414+054, B1030+074, 
J1632-0033, J1838-3427 and B1938+666, each of which 
is known to be strongly lensed by a single early-type galaxy, 
and B0128+437, B0739+366 and B1555+375 whose lens 



types are unknown. For the last three sources we take wj, = 
0.8 while Wk = 1 for the rest. Notice that the lensed sources 
of the CLASS statistical sample are not included in the like- 
lihood of the relative image separation probabilities because 
the likelihood of the CLASS statistical sample includes the 
relative image separation probabi lities as we ll as the abso- 
lute lensing probability (see, e.g., Chac 2003). Fig. [3] shows 
the confidence limits in the plane of /i (equation(6]l and v2/o 
(equation [10]). Here we have chosen ctsieo = 200 km s _1 for 
the fiducial velocity dispersion. Fig. [4] shows how the fac- 
tor y / 2/(= \/2o'siE/t>e,vir) behaves as the velocity dispersion 
(csie) is varied. Notice that v^csie would be equal to the 
circular velocity in the optical region if the galaxy mass pro- 
file were isothermal up to a few effective radii. The quantity 
v2csie will thus be our estimate of the circular velocity at 
the effective radius v c ,ett from strong lensing based on the 
SIE model (equation I12[l . To understand better Fig. [4] it is 
useful to rewrite / as 



f = fo 



CSIE0 



(16) 



which can be derived from equation §6§ using equation (|10[) . 
From equation (116[) we can see that for the given uncertain- 
ties of jj, and fo (Fig. [3} the uncertainty of / increases as 
asm decreases since fj, is negative. We can equivalently say 
that the uncertainty of the ratio / = usiE/'Uc,vir increases 
as trsiE (and thus u c ,vir) decreases for a given uncertainty of 
v c ,vii- Finally, Fig. [5] (Fig. [6j shows how f c , v ir (csie) behaves 
as (tsie {vc,vit) is varied. 

Next we consider the case where we use only the first 
part in equation (|lip . namely the relative image separa- 
tion probabilities. Notice that in this case the absolute 
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Figure 3. Likelihood contours in the plane of fi and \/2/o based 
on strong lensing and semi-analytical model of galaxy formation. 
Here parameter fi is the power-law index in the model f(v c v i r ) 
(equation [6]l for the ratio of the velocity dispersion (in the opti- 
cal region and the inner halo) to the virial circular velocity and 
parameter /o is the value of the ratio at the fiducial velocity dis- 
persion of ctsieo = 200 km s — 1 . The red contours are based only 
on the image separation distribution of radio lenses. The green 
contours are based on both the image separation distribution of 
radio lenses and all the lensing properties (including the lensing 
rate) of the CLASS statistical sample. The two contours for each 
case correspond respectively to the 68% and 95% confidence lim- 
its. 



abundances of the early-type galaxies from the SAM have 
no effects on the constraints on the velocity profiles of 
gal axies. For th is case we use all the lens systems used 
by I Chad ([2005), namely the 8 lens systems above along 
with the following 7 lens systems in the CLASS statistical 
sample; B0445+123, B0631+519, B0712+472, B1152+199, 
B1422+231, B1933+503, and B2319+051. The results for 
this case are also displayed in Figures [3] to [6] 

Bearing in mind that the strict range of validity for our 
analyses is 200km s^ 1 < « c . v ir J$ 600km s _1 (see above), 
our main findings are as follows. First, parameter [i must be 
negative meaning that the ratio of the inner velocity to the 
virial velocity must be larger for a less massive halo (Fig. [3]). 
A constant ratio (i.e. — 0) is excluded at a highly signifi- 
cant level. We have fi = -0.78±° ° l (68% CL) based on both 
the image separation distribution of radio lenses and the 
statistics of the CLASS statistical sample or ^ = — 0.67lg jj8 
(68% CL) based only on the image separation distribution 
of radio lenses. Second, the inner velocity dispersion ctsie 
and the virial circular velocity t> ClV ir scale as 



200 km s" 



200 km s" 



1.17_ 



u.^_ 0i04 



' 200 km s" 



for the case of using both the image separation distribution 
and CLASS statistics, or 



CTSIE 



200 km S" 



??CTccnt 

200 km s- 



0.82Z 



' 200 km s- 



,+o.n 

-0.08 



for the case of using only the image separation distribution. 
We take 



CTSIE 
& cent 



1.0±0.1 



(19) 



from IChad (2005) and lTreu fc Koopmansl ([2004) . Third, for 
galaxies with velocity dispersion ctsie 1^ 210-230 kins" 1 
the ratio v2o"siE/«c,vir(= Vc,eff/«e,vir) becomes increasingly 
larger than 1 as ctsie decreases (Fig. [4| . For a typical bright 
galaxy with ctsieo = 200 km s -1 , x/2/o = V^crsiEo/wc.virO = 
1.65^037 (image separation distribution + CLASS statis- 
tics) or 1.16lo 3o (i ma g e separation distribution only) at 
the 68% confidence level. However, for large galaxies with 
ctsie > 210-230 km s _1 , the ratio \/2CTsiE/i>c,vir < 1. For 
a galaxy with ctsie = 260 km s _1 , which corresponds ap- 
proximately to w e ,vir ~ 550km s _1 (see Fig. [5] and Fig. [6]), 



V2f = V2ct S ieA 



(image separation distri- 



bution + CLASS statistics) or 0.68jlo'ii (image separation 
distribution only) at the 68% confidence level. We discuss 
appropriate interpretations of these results along with pos- 
sible sources of systematic errors in §4. 



4 DISCUSSION 

In this work we have compared the velocity dispersion at 
about the effective radius of the optical region (ctsie) de- 
rived from strong lensing based on the SIE lens model (equa- 
tion [12]) with the circular velocity at the virial radius of the 
surrounding halo (w c ,vir) predicted by semi-analytical stud- 
ies of galaxy formation for the early-type galaxy popula- 
tion. Assuming that the inner velocity dispersion is related 
to the virial circular velocity by the simple power-law model 
(equation [6}, we find that the ratio of the velocity disper- 
sion to the circular velocity becomes increasingly larger as 
the velocity dispersion decreases; the power-law index n in 
equation ((6} must be negative (Fig. [3}. Hence, we can con- 
clude that there is clearly the trend that the smaller the 
surrounding halo is, the more enhanced the optical veloc- 
ity dispersion is. This is the most robust result from our 
work. In addition to thi s , lensing stu dies l|Kochaneklll99"i 
iTreu fc Koopmansl 12004 IChaej l2005t) have found that the 
velocity dispersion in the optical region derived from strong 
lensing based on singular isothermal ellipsoids, as is the case 
in this work, is about the same as the central stellar veloc- 
ity dispersion determined from spectroscopic observations, 
namely r\ = ct sie /ct cent ~ 1 (equation I19[) . 

Strong lensing gives the scaling between the central 
velocity dispersion CT cen t and the virial circular velocity 

22 +0.05 

fc,vir as CTcent oc v c vlr -0 ' 04 (image separation distribution 

0.33± 011 

+ CLASS statistics) or CT CO nt oc v c vir _0 08 (image separa- 
tion distribution only). Our results are broadly consistent 
with those from recent other studies b a sed on halo occu- 
pation statistics i|Vale fc Ostrikerl l2004 I Yang et al 1120051) 
and galaxy-galaxy weak lensing ( Mandelbaum et al. 120061 ). 
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Figure 4. The behaviour of v2/ (= V^csiE/^c.vir) as a function of crgiE- Here parameter ctsie is the velocity disp ersion (in the optical 
region and the inner halo) implied by strong lensing assuming the sing ular isothermal elli psoid model of[chac (2003) and parameter v c v i r 
is the circular velocity at the virial radius of the halo predicted by the lKang et alJ j2005T l semi-analytical model of galaxy formation. For 
our adopted model of the singular isothermal ellipsoid, the factor v^^SIE corresponds to the circular velocity in the optical region and 
the inner halo V ce g. The red line is based only on the image separation distribution of radio lenses. The green line is based on both the 
image separation distribution of radio lenses and all the lensing properties (including the lensing rate) of the CLASS statistical sample. 
The dashed and dotted lines represent respectively the 68% and 95% confidence limits. Notice that the uncertainty in the ratio becomes 
larger at smaller velocity for a given uncertainty of the velocity. 



These studies all obtained some relations between the virial 
mass of a halo M v - lr and the luminosity L c of a central 
galaxy hosted by the h alo. Suppose L c oc AfV 7 and the 
Faber- Jackson relation (jFaber fc Jackson] ll97q ) L c oc a 13 
where we shall take j3 = 4. It follows then <T ccnt oc ff^jf 
asMvir oc v 3 vir (see iBullock et alJlioOlh . I Vale fc Ostrikerl 
l|2004l ) found 7 ss 3.57 using halo occupation statistics for 
massive halos which implies a cont oc w^'vir- Using galaxy 
groups in the Tw o Degree Field Galaxy Redshift Survey, 
lYang et all (j2005T l found 7 w 4 for M vir > 10 13 /i _1 M Q but 
7 « 3/2 for M vir < 1O 13 /i _1 M . The range of the halo 
mass corresponding to our velocity range extends across 
the transition mass 10 13 M©, so the implied scaling will 
be somewhere between a ce nt oc v^' v f t and CT cen t oc Vr.'.v\r- 
The g alaxy- galaxy weak lensing study of lMandelbaum et al.l 
(2006) found 7 as 2.7 for galaxies with L c > L* which im- 
plies CTcont oc «e,vir- All these results are clearly consistent 
with our results within the statistical errors. 

Our plot (Fig. 0} of the ratio V / 2/(= \/2VsiE/«c,vir = 
w c ,ea/wc,vir) (equation [5] where ctsie is the velocity disper- 
sion at about the effective radius and is the circular- 
velocity at the virial radius) against the velocity disper- 
sion asm shows some interesting features. First of all, for 



galaxies with velocity dispersion ctsie < 210-230 km s _1 the 
ratio s/2f > 1. For example, a typical bright galaxy with 
ctsieo = 200 km s _1 has V%fo = l-65lo'.37 (image separa- 
tion distribution + CLASS statistics) or 1.161q'3q (i ma g e 
separation distribution only). Given v Cl<3 g = \/2a for our 
lens model (e quation 11211. the above results are in agree- 
ment with the lSeliakl (j2002h galaxy-galaxy weak lensing re- 
sult of fJ c ,opt/^c,vir = 1.68 ± 0.2 for an L* early-type galaxy 
with his adopted value of CT CGn t = 177 km s _1 |f| A striking 
feature of Fig. [4] is its behaviour at high velocity disper- 
sions. The ratio v2/ is about 1 or less for ctsie > 210- 
230 km s _1 . For a galaxy with ctsie = 260 km s^ 1 , which is 
hosted by a halo of v c ,vii ~ 550 km s _1 (see Fig. [6j|, our 
result is v2/ = 0.65^o '12 (image separation distribution 
+ CLASS statistics) or 0.68^°'^ (image separation distri- 
bution only) at 68% confidence level. The 95% confidence 
range is 0.42 ^ V2f ^ 0.99 (image separation distribution + 
CLASS statistics) or 0.48 ^ y/2f ^ 1.63 (image separation 



3 It appears that the value of V c> opt/v c ,viv changes little between 
Qcent = 177 km s — 1 and 200 km s _1 according to the results by 
ISeliakl <|2002t) . 
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Figure 5. The behaviour of function of ctsie- Here 

parameters v c y iT and cgiE are the same as in Fig. [4] The red and 
green lines are also the same as in Fig. [4] The dashed and dotted 
lines represent respectively the 68% and 95% confidence limits. 
The solid black line represents v c v ; r = \/2agiE. Notice that the 
uncertainty in v c v i r is nearly unchanged as <tsie is varied. 



distribution only) for asm = 260 km s" 1 . From lSeliakl (|2002l ) 
s/2f = 1.3 ± 0.2 for cr cent = 290 km s" 1 . Given that our de- 
rived velocity dispersion ctsie is believed to be similar to the 
central stella r veloci t y disp ersion cr ccnt fequation ll9[) . our re- 
sults and the lSeliakl 1)20021) results show some intriguing dif- 
ference although the 95% statistical errors partially overlap. 
For the large circular velocity systems (u c ,vir ~ 550 km s _1 ), 
it is likely that the haloes may host several galaxies and our 
inner velocity dispersion is appropriate for a central galaxy. 
One possible interpretation of our results is then that the 
optical galaxy is well within the peak radius r pea k where the 
circular velocity is highest and the baryonic boost of the 
velocity within the optical galaxy is relatively small for the 
halo with u c , v ir ~ 550 km s _1 (see below for further discus- 
sion). 

While our results appear to be in reaso nable agreement 
with halo occupation numb er studies of IVale fc Ostrikerj 
l|2004l and lYang et~aH d2005f) and the most recent weak lens- 
ing studies of lMandelbaum et al.l ([2006') , it is still important 
to consider possible sources of systematic errors for our re- 
sults. First of all, the abundances of early-type galaxies at 
low virial circular velocities from the SAM are uncertain at 
present. As can be seen in Fig. [TJ the abundances of early- 
type galaxies from the SAM appear to be abnormally high 
at low virial circular velocities a nd do not fit into our model 
CVF. In fact. lKang et al.l (120051 ) find that their predicted lu- 
minosity function of early-type galaxies does not match well 



observed luminosity functions at the low end of luminos- 
ity. Semi-analytical methods have their limitations primar- 
ily due to uncertainties in the star formation and feedback 
processes (which are also present in hydrodynamical simu- 
lations). In most semi-analytical studies it appears that the 
number of galaxies at the faint end tends to increase faster 
than the observed trend as the luminosity decreases. We 
have chosen to ignore the lowest four data points of the SAM 
(Table 2) in fitting our model CVF. As far as our analyses 
are limited to the range 200 km s _1 < w c , vir < 600 km s _1 , 
ignoring the lowest four data points would not affect the re- 
sults based on the image separation distribution only. How- 
ever, ignoring the lowest four data points may cause a sys- 
tematic error for the results based on both the image sepa- 
ration distribution and CLASS statistics because the lensing 
rate depends on the entire range of the velocity dispersion 
function but discounting the four data points may have bi- 
ased the behaviour of the circular velocity function so that 
the relation between the two may be affected. Nevertheless, 
it is likely that the systematic error might not be too large 
for the following reasons. Our procedure does not entirely 
ignore the low circular velocity galaxies but uses the extrap- 
olated abundances from the fitted function rather than the 
given abundances. Moreover, current observations show that 
as the luminosity of the early-type galaxy decreases below 
L» the surface brightness distribution becomes less concen- 
trated so that it is less effective for strong lensing. 

Secondly, the abundances of early-type galaxies at the 
highest velocity end might also be systematically biased. In 
fact, semi-analytical models tend to over-predict the number 
of very luminous galaxies, although the inclusion of AGN 
feedbacks seems t o cure this problem (|Bower et alj 120061 ; 
ICroton et al.|[2006l ). An over-predicted abundance of galax- 
ies at a given virial circular velocity will lead to an under- 
estimate of / at the corresponding velocity dispersion. To 
see this effect we have done the following numerical exper- 
iment. We decreased the numbers of early- type galaxies for 
bins i = 9 to 15 in Table 1 (with fJ ClV ir > 312.5 km s" 1 ) by 
successively larger proportions from i = 9 to 15 so that at 
i = 15 the adjusted number becomes one half of the unad- 
justed number. In this case we indeed find that the value of 
/ somewhat increases at large velocity dispersions in par- 
ticular at (Jsie = 260 km s _1 but within the 68% statistical 
error. 

Thirdly and finally, our results might be biased because 
of the assumed lens model, namely the SIE (equation I12p . 
Theoretically, the circular velocity (or the velocity disper- 
sion) derived from strong lensing depends on the radial pro- 
file and the shape of the model mass distribution. For the ob- 
served image size and morphology, the predicted size of the 
Einstein ring i?Ein and the projected mass within the ring 
MEin depend on the radial profile. However, for most of the 
known lens systems detailed mass modelling shows that the 
derived Rsin and MEin vary little as the radial profile is var- 
ied. The circular velocity at f?Ein depends on the mass within 
the sphere of radius i?Ein, Af(i?Em) rather than M^in- For a 
spherical mass distribution of the form p ~ r~ v a model with 
v — 1.7 (y = 2.3) would give a circular velocity at r = i?Ein 
12% lower (9% higher) than the isothermal model with v = 2 
for the same MEin- It is possible that the average mass pro- 
file of early-type galaxies systematically varies from smaller 
systems to larger systems. In fact, our own results appear to 
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Figure 6. The behaviour of <tsie as a function of i> CjV ir- Notations are the same as in Fig. [5] 



imply varying structures of early- type galaxies (see below). 
Nevertheless, the errors for / arising from approximating all 
the galaxies using isothermal models appear to be within the 
current 68% statistical errors. Furthermore, the inner mass 
profiles of early-type galaxies at least up to the effective 
radii appear to be close to isothermal from stellar dyna mi- 
cal modelling (e.g- lRix etafl 1 19971 ; [Thomas et al.ll2005l ) as 
well as from detailed modelling of individual lenses (e.g., 
iRusin fc Kochanekll2005l ; iKoopmans et al] |2006f). Thus, the 
assumption of the SIE model might not actually cause as 
large errors as estimated above. The shape of the mass dis- 
tribution also affects the derived circular velocity (or ve- 
locity dispersion) from strong lensing. For example, p rolate 
and o blate isothermal models give different results fsee lChael 
2003). In this work we have assumed that one half of galax- 
ies are oblate and the other half are prolate. The average 
shape of early-type galaxies might systematically vary from 
smaller s ystems to lar ger systems. However, according to the 
model of I Chad (|2003l ) the change of the velocity dispersion 
from the half oblate and half prolate case is ~ 10% even 
for the extreme cases of all oblate and all prolate. Consid- 
ering the above estimates of possible systematic errors, we 
can tentatively conclude that the essential trend of / for 
200km s _1 < u c , vir < 600km s~ x and the result \/2f < 1 
at ctsie = 260 km s _1 are likely to be real. 

The fact that \/2/(= « c ,eff/wo,vir) is not unity shows 
that galaxies are not well approximated by isothermal pro- 
files for significant parts of the virial radii. The trend of / as 
a function of the velocity dispersion constrained from strong 



lensing in this work is in agreement with theoretical expec- 
tations. F irst, more massive haloes have smaller co ncentra- 
tions (see iNavarro et al.l [l997l ; iBullock et all 120011 ). and so 
they have smaller peak velocities relative to the virial circu- 
lar velocity compared with smaller haloes. Second, the bary- 
onic modifications to the inner mass profiles of the haloes are 
less pronounced for more massive haloes and consequently 
the baryonic boost of the velocity may be smaller. On the 
other hand, smaller circular velocity systems may have a 
larger value of / as baryonic cooling could be more efficient 
in such systems (although feedback processes may suppress 
cooling in very small systems). Third, for massive haloes the 
radius of the optical galaxy may be much smaller than the 
peak radius r pca k of the halo where the circular velocity ob- 
tains its maximum while the optical radius is comparable to 
)"peak f° r small haloes. However, to assess the full significance 
of our results, detailed modelling is necessary incorporating 
theretical halo profiles, baryonic modifications of the haloes, 
observed light distributions and realistic d istribution func- 
tions (see, e.g., Ivan der Marel et~aT] l2000l for modelling of 
galaxy clusters). 

In this work, we consider only a si mple two-component 
mode l recently used in the literature l|Keetonll200l ; ISeliakl 
2002) to see whether our derived ratio of the SIE veloc- 
ity dispersion to the virial circular velocity can be matched 
qualitatively and shed new light on galactic structures. (We 
postpone a more detailed work to a future publication.) 
The model starts from an initial N FW halo and then turns 
into a Hernquist l|Hernquisj [l990') light (stellar mass) dis- 
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Figure 7. The circular velocity profiles for a two-component early-type galaxy model at z = 0.6 consisting of a Hcrnquist light (stellar 
mass) distribution residing in a halo, which is supposed to form from an initial NFW halo by baryonic cooling. The solid, dashed, 
dash-dotted, and dotted lines represent respectively the profiles for the final total mass, initial NFW halo, final modified halo, and 
Hernquist stellar mass distributions. For the lower panel (v cv i t = 200km s — 1 ), we take c v j r = 9 and F coo i = .09. For the upper panel 
(^c.vir = 550 km s _1 ), we take c v ; r = 5 and F coo i = 0.01. Here c v ; r is the halo concentration I Bullock et al 1l200lh and F coo i is the fraction 
of cooled baryons forming stars out of the initial total mass. In these models, the predicted ratios of the circular velocities at the effective 
and the virial radii v c e ff/v c v ; r are consistent with our derived values of ^2asis/v c V \ T . 



tribution by baryonic infall residing in a modified halo. We 
use th e adiabatic contraction model by iBlumenthal et all 
( 1986) to calculate the modified halo profile. Fig. [7] shows 
the constructed models at redshift z = 0.6 for the haloes 
of Uc.vir = 200km s _1 and « c ,vir = 550km s" 1 . For the 
v c ,vir = 200 km s _1 halo, we take c v i r = 9 (where cVir is 
the halo concentration) and -Fcod = 0.09 (where -F CO oi is 
the the cooled baryon fraction out of the total mass). For 
the n c ,vir = 550 km s -1 halo, c v ir = 5 and F coo i = 0.01. 
Here we have (approximately) estimated the values of c v i r 
using the simulation results of iBullock et all (|200ll ) while 
the values of F coo \ have been adjusted to produce circular 
velocity functions that are consistent with the results from 
this work. Consequently, in these models the predicted val- 
ues of the ratio of the circular velocities at the effective and 
the virial radii v c , c ff /i> CjV ir are in agreement with our derived 



values of y/2f for ctsie = 200 km s" 1 and 260 km s" 1 . The 
circular velocity curves shown in Fig. [7] imply evidently the 
following for galaxy formation processes and resultant mass 
profiles. First, the galaxies with u c , v ir = 200 km s _1 (and 
corresponding asm = 200 km s _1 ) and u CjV ir = 550 km s _1 
(o"sie = 260 km s _1 ) require very different values of Fcooi 
in order to be consistent with our derived values of y/2f. 
The fitted values of .F C ooi imply that the fraction of infalling 
baryons forming stars and consequential modification to the 
halo are becoming increasingly smaller as the halo gets more 
and more massive. This is qualitatively consistent with the 
fact that for galaxy clusters (the high mass limit of haloes) 
the NFW halo profile is well preserved. Second, the system of 
^ c ,vir = 200 km s _1 has a density profile close to isothermal 
up to a few effective radii but steeper than isothermal be- 
yond. On the other hand, the system of u c , v ir = 550 km s" 1 
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has a profile shallower than isothermal up to about a half 
of the virial radius. However, more detailed investigations 
along with more precise interpretations will be considered 
in a future publication. 

The parameterisation we use in equation ((SJl might not 
be optimal. The current CLASS and PANELS samples are, 
however, too small to allow us to explore a more realistic 
functional form. However, one advantage of the method is 
worth emphasizing: the lensing probability scales as / 4 and 
the separation scales as f 2 - the lensing properties therefore 
depends on quite high powers of /. This also implies that 
even if the semi-analytical modelling is somewhat uncertain, 
we may still be able to put strong limits on this parameter 
with the next-generation lens surveys. Current generation 
hydrodynamical simulations cannot yet resolve and simulate 
the inner parts of early-type galaxies realistically. For exam- 
ple J _the_liig2iresolut ion simulation of an early- type galaxy 
bv lMeza et all (120031 ) predicts a central velocity dispersion 
as high as 650 km s" 1 due to the compact size of their stellar 
component. So gravitational lensing can play an important 
role in empirically assessing the roles of baryonic cooling in 
galaxy formation and evolution. 
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